Level Spacings in Random Matrix Theory and Coulomb Blockade Peaks in Quantum Dots 



oo 

o 

O 

(N 

X> 
^1- 



i 

-a 
c 

o 
o 



> 

o 

(N 

o 
i> 
o 



X 



Damir Herman, T. Tzen Ong, 2 '0 Gonzalo Usaj, 2 ' 3 Harsh Mathur, 1 and Harold U. Baranger 2 

'Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106-7079 
2 Department of Physics, Duke University, Box 90305, Durham, North Carolina 27708 
Instituto Balseiro and Centra Atomico Bariloche, Comision Nacional de Energia Atomica, (8400) Bariloche, Argentina. 

(Dated: August 17, 2007) 

We obtain analytic formulae for the spacing between conductance peaks in the Coulomb blockade regime, 
based on the universal Hamiltonian model of quantum dots. New random matrix theory results are developed 
in order to treat correlations between two and three consecutive spacings in the energy level spectrum. These 
are generalizations of the Wigner surmise for the probability distribution of single level spacing. The analytic 
formulae are shown to be in good agreement with numerical evaluation. 
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I. INTRODUCTION 

The Coulomb blockade of electrons has been a remark- 
able tool for probing fundamental physics in nanoscale 



systems 



1.2.3.4 



In a semiconductor quantum dot (QD), weakly 



coupled to leads via tunneling, the effect is manifested as a 
series of spikes in the conductance of the device as a func- 
tion of a gate voltage which controls the number of electrons 
on the QD. In the quantum Coulomb blockade regime (de- 
fined below), the spacing between the spikes is determined by 
the ground state energy of the QD. Experimentally, the spac- 
ing shows mesoscopic fluctuations, and there have been sev- 
eral detailed studies of the statistical distribution of the peak 
spacings, 5 . 6 - 7 - 8 . 9 . 10 - 11 - 12 

The distribution of peak spacings results from the inter- 
play of electron-electron interaction and randomness. A typ- 
ical semiconductor QD in the experiments of Ref. 6 consists 
of a droplet of several hundred electrons confined to a two 
dimensional region of a few tenths of a micron in size. In 
the absence of electron-electron interaction, the motion of the 
electrons would be randomized by impurity scattering or by 
chaotic scattering from the boundaries of the QD; hence, it 
is expected 1314 that the single particle energy levels in the 
absence of interaction would be described by random matrix 
theory (RMT)^ In a real QD, however, it is necessary to take 
the additional effects of interactions into account. This is in- 
dicated by the failure of (essentially) non-interacting models 
to account for the observed distribution of Coulomb blockade 
peak spacings. 5 - 6 

The interplay of electron-electron interaction and random- 
ness is a notoriously hard problem. Nonetheless, due to the 
finite size of the QD, it has been possible to make significant 
progress. In particular, it is now believed that the energy lev- 
els of a weakly interacting QD are described statistically by 
a "universal Hamiltonian" (see Refs.d anrfl and references 
therein). According to this model, each QD is characterized 
by a set of single particle orbitals with single particle energies 
and wavefunctions distributed according to the appropriate en- 
semble of RMT. Given the occupation numbers of these single 
particle orbitals, the non-interacting contribution to the energy 
of a QD follows directly. The interaction contribution is de- 
termined entirely by the net charge on the QD (the "charging 



energy") and its total spin (the "exchange energy"). Thus the 
universal Hamiltonian is 



(1) 



where hi„ = or 1 is the occupation number of orbital i with 
spin a =t or J,, 5N is the number of excess electrons on 
the QD, and S is the total spin-operator. C denotes the ca- 
pacitance of the QD, and J is the exchange constant. Strictly, 
there is a fourth term in the universal Hamiltonian correspond- 
ing to interaction in the Cooper channel. However, this term 
vanishes when time reversal symmetry is broken by the ap- 
plication of a magnetic field, and, even when present, it is 
significantly smaller than the others and may be neglected. 3 

It should be emphasized that the universal Hamiltonian 
does not help calculate the energy levels of any particular QD. 
Rather it provides a description of the universal statistical fea- 
tures of the levels of all QDs that belong to the universality 
class under consideration. 

The form of the universal Hamiltonian is dictated by the 
few symmetries that remain to a random system such as a 
QD»i£ For example, Eq. (fl} applies provided the QD has spin- 
rotation invariance (no spin-orbit or spin-flip scattering). If 
in addition time reversal symmetry is intact, the single parti- 
cle levels are distributed according to the orthogonal ensemble 
of RMT; if time reversal symmetry is broken (by a magnetic 
field, for example), the unitary ensemble applies. 

If the exchange constant is set to zero, the universal Hamil- 
tonian reduces to the old constant interaction model which has 
proven unsuccessful in accounting for the observed peak spac- 
ing distribution. It is instructive to compare the physics of 
the two models. According to the constant interaction model, 
each filled QD level is doubly occupied in the ground state, 
except for the top level, which would be singly occupied if 
the total number of electrons in the dot were odd. Thus the 
total spin of the ground state is zero or one-half depending 
on whether the number of electrons is even or odd. Within 
the universal Hamiltonian model, higher spin ground states 
are possible: the higher single particle energy of these states 
is offset by the exchange energy which favors parallel align- 
ment of spins»i2ii& Indeed if the exchange constant is above a 
certain threshold, the QD should have a ferromagnetic ground 
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state that is fully spin-polarizedji&i^ Here we restrict our at- 
tention to the paramagnetic regime in which the ground state 
may have a substantial spin but is not fully polarized. 

It is also instructive to compare the physics of QDs to 
that of atoms. Atoms can be approximately understood by 
considering that the electrons occupy a set of self-consistent 
hydrogen-like orbitals. The precise filling of the incomplete 
shells (which determines the total spin of the atomic ground 
state) is then controlled by the interplay of electron-electron 
interaction and spin-orbit coupling encoded in Hund's rules. 
The universal Hamiltonian provides a comparable description 
of the energy levels of weakly interacting QDs. The key dif- 
ference between atoms and QDs is that for atoms the confin- 
ing potential is spherically symmetric whereas for QDs it is 
essentially random. 

The universal Hamiltonian is found to give a good quan- 
titative account of the experiments of Ref. 6 when finite 
temperature effects and small non-universal corrections are 
included . 1 2 ' 20 ' 2 1 - 22 i 23 Here, we obtain analytic formulae for the 
distribution of conductance peak spacings in the low tem- 
perature and large QD limit. Though the distribution has 
been evaluated numerically before ^ 21 ' 22 ' 23 analytic expres- 
sions should prove helpful in the analysis of experimental 
data. 

The quest to develop such analytic expressions raises some 
interesting problems in RMT Within the constant interaction 
model, the distribution of peak spacings, II(s), is controlled 
by the distribution of the spacing between consecutive energy 
levels, denoted Pi (A). Level spacing distributions are noto- 
riously hard to calculate but by now this distribution is well 
understood. 15 The exact analytic expression for -Pi (A) is far 
too complicated to be useful in practice. Fortunately, the cel- 
ebrated Wigner surmise provides a simple and remarkably ac- 
curate approximation to the exact result. Within the universal 
Hamiltonian model, however, the peak spacing distribution is 
controlled by the joint probability distribution of the spacings 
between several consecutive levels, which we shall denote as 
P 2 (Ai, A 2 ), P 3 (Ai, A 2 , A3), etc. As for the single spacing 
distribution, the exact expressions for the joint spacing distri- 
butions are too cumbersome for practical use. 

We have, therefore, developed approximate expressions 
analogous to Wigner's surmise for the joint level spacing dis- 
tributions and used these to obtain formulae for the conduc- 
tance peak spacing distribution. In addition, we have de- 
veloped new numerical techniques for evaluating the exact 
expressions which are much more efficient than previously 
known methods. The approximate analytic expressions are 
in excellent agreement with the numerics. 

To our knowledge there has been no previous investiga- 
tion of these joint spacing distributions in a physics context. 24 
Number theorists have studied a quantity dubbed the nearest 
neighbor spacing distribution in connection with the zeros of 
Riemann's zeta function. 25 That distribution is closely related 
to P 2 (Ai, A 2 ), and, indeed, the zeta function work, as well as 
other known results, provide a useful test of the accuracy of 
the approximate formulae we develop. 

The rest of the paper is organized as follows. First, we start 
with the main results of interest from a QD point of view: 



Section II gives our analytic results for the peak spacing dis- 
tribution. Then, in Section III we discuss the RMT quantities 
on which the peak spacing distribution is based. The approx- 
imate expressions for the joint distribution of several consec- 
utive levels are developed here. Section IV presents our exact 
results for these RMT quantities and compares them with the 
previously given surmises. Finally, in Section V we summa- 
rize and conclude. 



II. PEAK SPACING DISTRIBUTION 

The peaks in the conductance occur at values of the gate 
voltage for which the number of electrons on the dot changes 
by one. In the quantum Coulomb blockade regime, the spac- 
ing between the N — 1 — > N and N — > N + 1 peaks is 

82 = Pat+i+Pat-1 — 2Pjy , (2) 

where En denotes the ground state energy of the QD with N 
electrons. The sequence of peak spacings is sometimes called 
the addition spectrum of the QD. For later use it is convenient 
to define the shifted peak spacing 

e 2 

8 = 5*- — . (3) 

It is instructive to first consider the constant interaction 
model. For N even, the N— 1 — > N transition involves adding 
an electron to level N/2, while in the transition an 

electron is added to level N/2 + 1. Thus, for an odd-even-odd 
transition (for brevity, an even spacing hereafter) 

s = A, (4) 

where A = e^r / 2 +i ~£n/2 is the spacing between the two lev- 
els. Similarly, for an even-odd-even transition (an odd spacing 
hereafter) 

s = 0. (5) 

Thus in the constant interaction model, the peak spacing dis- 
tribution is bimodal: the even spacings are distributed in the 
same way as the single particle level spacings while the odd 
spacing distribution is a delta function spike at s = 0. 

Let us now consider the case of the universal Hamiltonian 
model, Eq. (fl}. For even N, the state with the lowest N/2 
single particle levels doubly occupied has total spin Sjv = 0. 
It is therefore an eigenstate of the universal Hamiltonian al- 
though, as we shall see, it is not necessarily the ground state. 
Next, consider the four states obtained by promoting an elec- 
tron from the highest occupied level to the lowest unoccupied 
one (the occupations here are stated relative to the constant 
interaction ground state). The single-particle energy of these 
states is greater by A, and they can be combined into a spin 
triplet and singlet. We have thus identified five eigenstates 
of the universal Hamiltonian: the constant interaction ground 
state with energy Eq, the degenerate triplet states with energy 
Eq + A — 2 J and the singlet state with energy Eq + A. Pro- 
ceeding in this way we can construct all the eigenstates of the 
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universal Hamiltonian from the excitations of the constant in- 
teraction model. 

For J = 0, the constant interaction ground state is the true 
ground state, but for J > A/2 the triplet states become lower 
in energy. In principle a state with still higher spin might have 
even lower energy: such states are costlier in terms of sin- 
gle particle energy but have large (negative) exchange contri- 
butions to their energy. One might expect that in the para- 
magnetic limit J < (A) (the mean level spacing) the single 
particle energy would dominate, and a high spin state would 
be unlikely. Indeed for J/ (A) as large as 0.5 we have veri- 
fied that, in systems without time-reversal symmetry, approx- 
imately 98% of the ground states for even N have Sjf = or 
Sn = 1. Similarly, Sn = \ and Sn = § are the dominant spin 
states for N odd. Hence, in the following we will consider 
only the competition between the two lowest spin states. 

Now consider the peak spacing for an even spacing. For 
simplicity, first assume that Sjv— l = Sn+i = \ and Sn = or 
1; in this case the spacing is given by 



A-§J forA>2J 



§J-A forA<2J. 



(6) 



Note that s > J/2. The cumulative probability that the spacing 
is smaller than a certain value is 







for s< ~ J 



\J+S 



F{s)= P l( A) d A for §J< S <§J (?) 



iJ+s 



Pi(A)dA for \J<8. 



Here Pi (A) denotes the single particle level spacing distribu- 
tion, for which both the exact result and an extremely accurate 
approximation — Wigner's surmise — are known in RMT^£ By 
differentiation we can convert the cumulative distribution to 
the peak spacing probability distribution 



>(*) = 



o 



Pi (U- 



-P (U- 



for s< i J 



for ±J<s< §J 



for §J<s 



Pi (*+§J) 

(8) 

Note that there is a sharp jump in the distribution n ovon at the 
left edge of its support (s = J/2), in contrast to the smooth 
behavior in the constant interaction model. Also n cvcn is con- 
tinuous at 5 J/2 but kinked in systems with time-reversal sym- 
metry [see Eq. (TT2l below]. 

The odd spacing is similar but more tedious to analyze be- 
cause there are four cases to consider: the QD may have either 
spin zero or spin one in both initial and final states. We de- 
note the spacing between the [(N— l)/2] th level and the level 
above it Ai and between the [(N + l)/2] th level and the level 



^ / 



above that A2. Straightforward analysis then shows that 
for Ai>2Jand A 2 >2J 

for Ai>2Jand A 2 <2J 

for Ai < 2 J and A 2 > 2 J 



s = < 



A 2 -±J 
Ax-ij 



(9) 



Ai + A 2 -|j for Ai<2Jand A 2 <2J 



Note -§J<s<§J. Provided Ai > 2 J and A 2 > 2 J (always 
true in the constant interaction limit of J — * 0) all odd peak 
spacings have the same shifted value s = § J independent of 
the values of the level spacings. However if these conditions 
are not met the spacing does depend on the values of Ai and 
A 2 . A short calculation reveals 

n oAd (s)=S{s--j) / P 2 (Ai,A 2 )tfAidA 2 +C7( s ). 

(10) 

Here P 2 (Ai, A 2 ) is the joint probability distribution of two 
consecutive level spacings, for which we will derive an accu- 
rate approximation analogous to the Wigner surmise in Sec- 
tion |nl] The continuum distribution is given by 



r 



for s<-|j or s>|j 



C(s)={ 



P 2 ([M, p) dv for-|j<s<-ij 



J-s 



P 2 (fx, ft) dv+ 2 J P 2 [v,s+- ) dv 



J 



for 



ij<s<§J 



(11) 

with fi=jj+j(s+v) and Jx = \j+\{s—v). 

Eqs. dH), ( I lot and (fTTT i are the expressions for the Coulomb 
blockade peak spacing distributions when it is assumed that 
the QD cannot have a ground state spin greater than one. 
These expressions apply whether time-reversal symmetry is 
intact or broken: the only difference in the two cases is that 
it is necessary to use the level spacing distributions Pi (A), 
P 2 (Ai,A 2 ) for the Gaussian orthogonal ensemble (GOE) 
when time-reversal symmetry is intact and for the Gaussian 
unitary ensemble (GUE) when it is broken. 

Unfortunately, this very simple approach is not enough to 
describe some important features of the peak spacing distribu- 
tion for values of J/ (A) ~ 0.3. States with S = | have to be 
included. Table U shows all possible spin transitions, the cor- 
responding peak spacing, and the conditions on the level spac- 
ings for the states involved to be the ground states. We used 
the spin selection rule AS'=±i upon the addition of an elec- 
tron. By a simple inspection of Table H] we note several new 
features that appear in this analysis: (i) The even peak spacing 
distribution now depends not only on Pi (A) and P 2 (Ai , A 2 ) 
but also on the joint probability distribution of three consec- 
utive level spacings, P3(Ai,A 2 , A3) — an accurate approxi- 
mation to this quantity is derived in Section [Till (ii) even tran- 
sitions with s < h J and odd transitions with s > | J are now 
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(Sn-1, Sn, Sn+i) 


Spacing 


Condition 


(1.0, i) 


s = A 2 -§J 


Ai+A 2 >3J 
A 2 >2J 




|J< S 


A 2 +A 3 >3J 




s = -A 2 + §J 


Ai+A 2 >3J 
A 2 <2J 




ij<s<|j 


A 2 +A 3 >3J 


(- 1 


s = A 3 -iJ 


Ai+A 2 >3J 
A 2 <2J 




-\J<s<lJ 


A 2 +A 3 <3J 


(- 1 -) 

V 2 ' ' 2> 


s = Ai-ij 


Ai+A 2 <3J 
A 2 <2J 




-|J< S <|J 


A 2 +A 3 >3J 


(- 1 -) 

V 2 ' ' 2 / 


s = Ai + A 2 +A 3 -|J 


Ai+A 2 <3J 
A 2 <2J 




-|J< S <|J 


A 2 +A 3 <3J 








(0,1,0) 


S =|J 


Ai>2J 
A 2 >2J 
Ai+A 2 >3J 


(o,|,i) 


s=A 2 -|J 


Ai>2J 
A 2 <2J 






Ai+A 2 >3J 


(1,1,0) 


s = Ai-|J 


Ai<2J 
A 2 >2J 






Ai+A 2 >3J 


(1,1,1) 


s = Ai+A 2 -§J 


Ai<2J 
A 2 <2J 




ij<s<|j 


Ai+A 2 >3J 


(l,f,l) 


s = -Ai-A 2 + |J 


Ai<2J 
A 2 <2J 




1J< S <|J 


Ai+A 2 <3J 



Even 



Odd 



TABLE I: List of all possible ground state spin-transitions up to 
S =% and their corresponding peak spacing. Only transitions with 
AS = ±| are considered. Last column shows the conditions on the 
nearest level spacings so that the states involved in the transition are 
the ground states. Note that the even cases involve three consecutive 
level spacings. 
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FIG. 1: Partial contributions to the peak spacing distribution calcu- 
lated using expressions in Appendix [A] for the GOE and for J = 
0.325(A). The distributions in the lower panels were multiplied by 
the indicated factor for the sake of comparison. Note the transitions 
involving 5* = | add tails to both the even and the odd distribution 
and add a discontinuity to the latter at s — \ J. 



involved the diagonalization of 10 6 random matrices of size 
100 x 100. Each random matrix corresponds to a different re- 
alization of the single particle levels in Eq. (1), and the peak 
spacing distribution is inferred by histogramming the calcu- 
lated peak spacings for these "virtual quantum dots". As can 
be seen from the figures, the agreement between the simula- 
tion and the analytic formula is excellent. 

We have found that it is not necessary to consider higher 
spins for J/(A) as large as 0.4. The analytic expression of- 
fers considerable advantage over the numerical calculation if, 
for example, it is necessary to compare the universal Hamilto- 
nian model to an experimental distribution with the exchange 
constant as fitting parameter. 



III. LEVEL SPACING DISTRIBUTIONS: SIMPLE 
APPROXIMATIONS 



possible; and (iii) a discontinuity appears in the odd distribu- 
tion at s= i J. 

An explicit expression for the peak spacing distribution is 
also possible in this case. We present it in Appendix lAl as the 
calculation is rather lengthy. Here, we simply show the results 
for the GOE in Figure Q] Transitions (|,l,§)and(§,l,|) 
both give the same contribution, so only one is displayed; sim- 
ilarly for (0, i, 1) and (1, ±, 0). 

A comparison of the peak spacing distribution, calculated 
in the approximation that allows no spin state higher than |, 
to a numerical simulation 20 is shown in Figures [2] and [3] for 
the GOE and GUE, respectively. The numerical calculation 



In this Section we develop accurate approximations to the 
joint probability density of consecutive level spacings. These 
expressions are needed to complete the expressions for the 
Coulomb blockade peak spacing distribution discussed above. 

First it is helpful to recall Wigner's work. 15 We are inter- 
ested in Pi (A), the spacing between two consecutive levels 
near the center of the spectrum of a large M x M random 
matrix. The exact result, derived by a lengthy and intricate 
calculation, 15 involves an infinite product of eigenvalues of 
prolate spheroidal functions which are difficult to evaluate nu- 
merically. It is therefore of limited practical value. 

Now, following Wigner, consider a 2 x 2 real symmetric 
Gaussian random matrix. Such a matrix has two eigenvalues, 
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FIG. 2: (Color online) Comparison between the numerical result 
(dots) and the analytical result (solid line) for the peak spacing dis- 
tribution in the Gaussian orthogonal ensemble with J = 0.325(A). 
The upper (lower) panel corresponds to the even (odd) case. The 
agreement is remarkably good. 



and an elementary calculation shows that the spacing between 
them is distributed according to 

P^ s (A)^Acx P (-Ja 2 ) . (12) 

The scale is chosen so that the mean level spacing is unity, 

(A>=/ APf^ s (A)dA = l. (13) 
Jo 

Eq. ( fT2b is the celebrated Wigner surmise. For the GUE, the 
surmise is similarly derived by considering a 2 x 2 complex 
hermitian Gaussian random matrix, yielding 



n S ws(A) = §A 2 ex P f-lA 2 



(14) 



The Wigner surmise is known to be an excellent approxima- 
tion to the level spacing distribution for the orthogonal ensem- 
ble where it has been checked against numerical calculations 
and rigorous bounds (see page 157 of Ref.ll5l). For the uni- 
tary case, to our knowledge, the surmise has not been subject 
to comparably rigorous scrutiny, but it is generally believed 
to be very accurate 26 (see also Fig. 8 below where the unitary 
Wigner surmise is compared to our numerical calculation of 
the level spacing distribution). 

Let us turn to the joint probability distribution of two con- 
secutive level spacings, P2(Ai, A2). The exact result may 
be derived using the theory of Fredholm or Toeplitz deter- 
minants, as in Section IIVDI below. However, a more useful 
approximation may be derived by considering, in the spirit of 
the Wigner surmise, a 3 x 3 real symmetric Gaussian random 
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FIG. 3: (Color online) Same as in Fig.[2]but for the Gaussian unitary 
ensemble. The agreement here is even better due to the fact that the 
probability of large spin is lower in this case. 



matrix (for the GOE). Using the standard joint probability dis- 
tribution of the levels of a random matrix, we obtain 



P 2 s . 3 e X 3(A 1 ,A 2 ) = 4 X /— a^A 1 A 2 (A 1 + A 2 ) (15) 



x exp 



2a 



[A 2 +A 2 +AxA 2 ] 



with a = 27/ (8tt). The distribution has been normalized and 
the energy scale chosen so that (Ai) = (A 2 ) = 1. We con- 
jecture that Eq. ( [Tol l is an accurate approximation to the true 
joint spacing distribution for large GOE matrices. 

A similar result for the GUE may be derived by considering 
3x3 complex hermitian Gaussian random matrices. We obtain 



pr3 c x3 (A l7 A 2 ) = 



Ab 4 



x exp 



A 2 A 2 (A 1+ A 2 ) 2 
^[A 2 +A 2 + AiA 2 ] 



(16) 



with b = 729/(1287r), which we conjecture is a good approxi- 
mation to the true joint spacing distribution for the GUE. 

There is no a priori justification for either Wigner's surmise 
or for our analogous conjectures. To test our conjectures we 
subjected them to many checks, some of which we describe 
here. For definiteness we focus here on the unitary ensemble; 
but we have done similar tests in the orthogonal case. 

First, we note the relationship 

/>oo 

P 1 sue (A 1 )=/ Pf uc (A 1; A 2 )dA 2 (17) 
Jo 

between the exact single level spacing distribution Py and the 
exact joint spacing distribution P 2 . Performing the indicated 
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FIG. 4: (Color online). Comparison between the GUE Wigner sur- 
mise and the 3x3 approximation Eq. J 1 8t . Inset: Difference between 
the latter and the former. The agreement is excellent. 
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FIG. 5: (Color online) The nearest neighbor spacing distribution 
Pnn(A) calculated using the Toeplitz representation Eq. j61 | i with 
n = 200 (black points) is in good agreement with the generalized 
Wigner surmise approximation, Eq. j_ 1 b (blue curve). 



integral using our surmise Eq. (1161 . we obtain 

/>oo 

A g ;3 e x 3 (Ai) = / ^(Ai.AaJdAa 
Jo 



(18) 



where 



h{x) = -x 3 e-* 2 +-xe~ x2 +^Q-x 2 +x 4 ^)Eric{x) (19) 
I 



and Erfc(x) is the complementary error function. 

Fig.|4]shows a comparison between the Wigner surmise and 
the 3x3 approximation Eq. (fT8l ). The agreement is excellent 
over the entire region where the distribution has substantial 
weight and also for small spacings. In the tails, the relative 
error is larger but, of course, very small on an absolute scale. 
To be precise 99% of the weight of Pf^ s (A) lies in the range 
< A < 2.11. Over this range, Pf^s and P i S ; 3x3 disa S r ee by 
less than 3%. 

Another quantity related to P 2 (Ai,A 2 ) is the nearest 
neighbor spacing introduced by Forrester and Odlyzko in con- 
nection with studies of the zeros of the zeta function:— the 
nearest neighbor is the closer of the level just above and the 
level just below a given level. The distribution of the nearest 
neighbor spacing is, therefore, 



ps n uc (A) = / dA x ,5(A-A 2 )P| uc (A 1 ,A 2 )dA 2 
Jo Jo 

+ dAi / <S(A — Ai) Pf ue (Ai, A 2 ) dA 2 . 

JO J A t 



(20) 



If we substitute the 3x3 approximation ( fl6l l into Eq. d20b . 
we obtain the 3x3 approximation to the nearest neighbor 

62 g nn ((&/6)2 A) where 



distribution,^, 3 (A) 



9nn(x) 



27 



■x 2 e- 12x 



[l8x + 84x 3 



+e 9a;2 7r2(3_4 a ; 2 + 4 x 4 )Erfc(3a;)] (21) 



Fig. [5] shows a comparison between this 3x3 approximation 
to P nn and the exact numerical computation in Section [TV] 
Again the agreement is very good except in the tails which 
have negligible weight. We have also confirmed that these two 
results are in good agreement with the exact numerical com- 
putation of P nn using the method of Forrester and Odlyzko. 25 

The covariance of consecutive spacings is also known.— 
The exact result for cov(Ai, A 2 ) is —0.922 in the orthogo- 
nal ensemble and —0.944 in the unitary ensemble; within our 
approximation for P 2 they are found to be — 87r/27~ —0.930 
and -32 (27 V3 - 8vr) /729 ~ -0.950, respectively. 

As a final check we can compare Eq. ( fl6] l directly to the 
exact numerical computation of P 2 discussed in Section [TV] 
Again the agreement is reasonable (see Fig.|9]below). 

Next, we turn to the joint probability distribution of three 
consecutive spacings. To this end, we consider a 4x4 real sym- 
metric Gaussian random matrix. From the standard expres- 
sion for the distribution of eigenvalues for a random matrix, 15 
we obtain the normalized distribution 



Pf^A^A^As 



=A 1 A 2 A 3 (A 1 +A 2 )(A 2 +A3)(A 1 +A 2 +A3)exp[--{2(A 1 +A 2 ) 2 +2(A 2 +A3) 2 ^A 1 +A 3 ) 2 }]. (22) 
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FIG. 6: (Color online). Comparison between the GOE Wigner sur- 
mise and the distribution obtained from the 4x4 approximation 
Eq. ( 1231 ) by integrating out either Ai and A2 (solid line) or Ai and 
A3 (dashed line). Inset: Difference between each approximation and 
the Wigner surmise. 



Using this distribution, we find the average of the middle spac- 
I 



ing, / = (A2) w 0.8388. A difficulty that now arises is 
that due to end effects in our extremely finite sized matrix, 
/' = (Ai) = (A 3 ) « 0.9400. We resolve this problem by the 
standard procedure of "unfolding"^ introducing the rescaled 
distribution 

P| 4 c x4 (A 1 ,A 2 ,A 3 ) = /' 2 /n 4 C x4(/'Ai,/A 2 ,/'A 3 ). 

(23) 

We conjecture that Pl^^Ai, A 2 , A3) is an accurate ap- 
proximation to the true joint spacing distribution for the GOE. 

We have subjected this conjecture to tests similar to those 
applied to the distribution P 2: 3x3 above. For example, we 
can integrate out the end spacings A x and A3 to obtain an 
approximation to the single level spacing distribution or we 
can integrate out the lower ones Ai and A 2 to obtain a second 
approximation. The two approximations are compared to each 
other and to the nearly exact Wigner surmise in Fig.|6]and are 
found to be accurate outside the tails. On the basis of tests 
such as these we conclude that our conjectured distribution 
Eq. d23l ) is sufficiently accurate for our purposes. 

The corresponding distribution for the unitary ensemble can 
be similarly derived by consideration of a 4 x 4 complex her- 
mitian matrix. We obtain the normalized distribution 



P| ; u 4 x4 (A 1 ,A 2 ,A3) = K A?A2A2(A 1 +A 2 ) 2 (A 2 +A3) 2 (A 1 +A 2 +A3) 2 exp[--{2(A 1 + A 2 ) s 



r 



-2(A 2 + A 3 ) 2 +(A 1+ A3) 2 }], 

m 



where k m 0.4789. In this case we find h' = (Ai) = (A3) w 
1.2288 whereas h = (A 2 ) w 1.1177. Unfolding the distribu- 
tion d24i l. we obtain 



P 3 s ; u 4 x 4 (Ai,A 2 ,A3) = /i' ; 



hA 2 ,h'A 3 ) 



(25) 

which we find to be an accurate approximation to the true joint 
spacing distribution for the GUE. 

In summary, we have obtained simple expressions for 
the joint spacing distribution of two consecutive spacings, 
Eqs. ( fT~6b and ( fToT i. and for three consecutive spacings, 
Eqs. d23l and d25l l. for the orthogonal and unitary ensembles, 
respectively. Together with the Wigner surmise, Eqs. (fT2b and 
O, these distributions are the random matrix theory tools 
needed to compute the Coulomb blockade peak spacing dis- 
tribution in the paramagnetic limit. Using these tools, we ob- 
tained the peak spacing distributions already given in Section 



IV. LEVEL SPACING DISTRIBUTIONS : EXACT RESULTS 

In this Section we derive a number of exact results in ran- 
dom matrix theory. These results allow us to certify the ac- 
curacy of the more useful simple approximations developed 
in the previous section. Some of the results, notably the new 



exact calculation of P nn (A) and the analysis of its asymptotic 
behavior, are also of intrinsic interest in random matrix theory. 

The problem under consideration is the following: We are 
given a set of n energy levels and P{x\, . . . , x n ), their joint 
probability distribution. According to random matrix theory 
the energy level distribution is 



P(xi 



-m n 

l<»<j'<n 



Xi 



tj I' 3 exp 



i=l 



(26) 

where the exponent j3 = 1, 2, or 4 corresponds respectively to 
the orthogonal, unitary, and symplectic ensembles of random 
matrix theory. M is a normalization constant, and 1 / is an 
energy scale set by the mean spacing between levels. 

^From the joint probability distribution P(x\, . . . , x n ) we 
wish to calculate quantities such as the consecutive level spac- 
ing distribution 

P\{t) = )Cin(n— 1) / dxz . . . / dx n 

J out J out 

P(x% — > --,x 2 ->■ -,x 3 ,...,x n ) , (27) 



where 



/—t/2 />oo 
dx + dx 
-00 Jt/2 



(28) 
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This is the probability that one level is at — t/2, another at 
t/2, and none of the others are in between. We may place 
any of the n levels at —t/2 and any of the remaining (n — 1) 
levels at t/2; this is the origin of the combinatoric pref actor 
in d27l >. The normalization K.\ ensures that J™dtPi(t) = 1. 
Generally we wish to calculate the level spacing distribution 
in the limit n — > oo. In addition to the level spacing, similar 
quantities of interest here are a more primitive quantity that we 
shall call the spacing determinant (defined below), the nearest 
neighbor spacing, and the joint probability distribution of two 
consecutive spacings. 

The essential difficulty in random matrix calculations is 
the highly correlated nature of the probability distribution 
Eq. d26j >. Level spacing calculations are particularly difficult 
due to the piecewise continuous nature of the integration do- 
main. Nonetheless, exact formal expressions for these quan- 
tities have been derived. For example, the spacing determi- 
nant can be expressed as the Fredholm determinant of a par- 
ticular integral operator K, the n — > oo limit of a certain 
n x n Toeplitz determinant, and the solution to an ordinary 
non-linear second-order differential equation. Explicit evalu- 
ation of the spacing determinant is then traditionally carried 
out numerically from the Fredholm determinant or the non- 
linear differential equation. In Section ITV Al we find that di- 
rect evaluation of the Toeplitz determinant for large but finite 
n provides a third efficient numerical method of calculating 
the spacing determinant. 

In Sections llV BIIIV Cl and llV Dl we derive new expressions 
for the consecutive level spacing Pi (A), the nearest neighbor 
spacing P nn (A), and the joint distribution of two consecu- 
tive spacings P 2 (Ai, A2), expressing each of these quantities 
as a Toeplitz determinant. These quantities may then be nu- 
merically evaluated by direct calculation of the determinants. 
By contrast, in the conventional approach they are expressed 
in terms of the eigenvalues and eigenfunctions of the kernel 
K, which are harder to handle numerically. The Toeplitz rep- 
resentation also makes it easy to obtain the large A asymp- 
totic behavior: the previously unknown asymptotic behavior 
of P nn (A) is given in Section ITV El For simplicity, we con- 
centrate upon the unitary ensemble, (3 = 2; the other ensem- 
bles will be studied in future work. 

Toeplitz determinants arise in many contexts, among them 
signal processing 27 and statistical mechanics. 28 Consequently, 
a great deal is known or conjectured about their asymptotic 
behavior and efficient numerical methods exist for their eval- 
uation; that is the chief virtue of the Toeplitz expressions for 
Pi (A), P„„(A), and P 2 (Ai, A 2 ) derived here. 



A. Level Spacing Determinant 

The level spacing determinant is the probability that a band 
of width t is entirely void of energy levels, 

E(t)= I dx x ... I dx n P(x u ...,x n ) , (29) 

J out J out 



where J out dx is defined in Eq. d28l ). Evidently E(0) = 1, 
whereas £'(oo)=0. By straightforward differentiation 29 



dE 



(30) 



', / dx 2 ... dx n P(xi — > — ~, x 2 , ■ ■ ■ , : 

J out •/ out ^ 



F is the probability that one level is at the edge of the band 
while the others lies outside it. Evidently, P(oo) = while 
P(0) = Pi(0) where Pi is the mean density of levels, or the 
one point correlation, defined as 



Ri(x) 



dx-). 



dx n P(xi 



x, x 2 ■ 



■ x n ^j 



(31) 

Differentiating once again, we find that the consecutive 
level spacing distribution, Eq. d27l i. is given by 



Pi(t)=K 1 



d 2 E 
dP 



(32) 



In deriving this equation we have assumed that the joint prob- 
ability distribution P(xi, . . . ,x n ) vanishes if two levels coin- 
cide, and have also neglected a bulk term in which the deriva- 
tive acts on the integrand in ( f3TT >. The error is found to vanish 
as n — ► 00. Eq. d32b shows that in principle Pi is determined 
by E. In practice, E is usually computed numerically, and 
hence Eq. ( f32l > is not the best way to determine Pi. How- 
ever, we may use this identity to show that the normalization 
constant is JC\ =1/Pi(0). 

It is now useful to briefly recount how the exact formal ex- 
pressions for E(t) are obtained. It is more convenient for this 
purpose to use Dyson's circular unitary ensemble rather than 
the Gaussian unitary ensemble [Eq. (126b with /3 = 2]. The cir- 
cular unitary ensemble describes n angles distributed around 
the unit circle according to the normalized distribution 



P(0i 



1 1 

n! (2tt)™ i 



n 



(33) 



<r<s<n 



It is well known that (up to an irrelevant scale factor that deter- 
mines the mean level spacing) the local statistical properties 
of these angles are identical to those of energy levels governed 
by the GUE. 

In the circular ensemble, E(t) is the probability that all an- 
gles lie outside the arc —t/2 < 9 < t/2. To compute E(t) it 
is helpful to rewrite 



P(0i,...,0») = ^£(-l) P (-l) c 



P,Q 



Xy>P(l)(0l)VQ(l)(0l) ■ ■ ■ <PP(n) {6n)V*Q(n) ( 6 n) 



(34) 



Here pi(0) = 1/v^tt, y> 2 (0) = e^/V^, . . . , y n {6) = 
e l ( n ~^ e /V2tt. P and Q represent permutations of the inte- 
gers {1, . . . , n} and (— l) p ^ is the parity of the permutation. 
This representation follows from (l33l by means of the Van- 
dermonde identity 



n 

l<i<j<n 



= det M. 



(35) 
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where M is the n x n matrix with elements M^, 
If we define 



v-i 



de Vr {e) V * s {6) = / exp[i(r- S )0], 

/t/2 Jt/ 2 

(36) 

it follows 



1 



' PQ 

—\ ^2i- l ) R 9P(l)MPW] ■ ■ ■ 9P(n),R[P(n)] 



PR 



5l,-R(l) • • • 9n,R(n)- 



(37) 



PR 



To obtain the second line, we have introduced Q = RP and 
used (— 1)Q = (— l) p (— l) fi . The third line follows from a 
simple rearrangement of the terms in the summand. The P 
sum is now trivial, and we find 



(38) 



E(t) = det g 



where g is an n x n matrix with matrix elements given by ( 

The mean spacing between angles is 27r/n (show by com- 
puting the one point correlation Ri). Hence we introduce A, 



(39) 



as a measure of the band-width in units of the mean level spac- 
ing. Making this change of variables, we may finally write 

E(A) = lim det g (40) 

n — >oo 

with g given by Eqs. d36b and d39l ). 
An n x n Toeplitz matrix has elements 

T rs = [ ^-f(0)exp[i(r-s)0] (41) 
Jo ^ 

that are controlled by a single function f(6). Comparing 
Eqs. ( l36l l and d4"TV we conclude that the matrix g is of the 
Toeplitz form with f(9) = l for irA/n < 9 < 2tt — nA/n and 
zero otherwise. Hence Eq. (l40l > is the representation of the 
level spacing determinant in terms of a Toeplitz determinant. 

To motivate the Fredholm representation, consider the 
eigenvectors of g, 



J29rs4 a) =X (a) 4 a) • (42) 

s=l 

Here a — 1, ... ,n labels distinct eigenvectors and their 
eigenvalues A^ Q ^. In the n^oo limit, we set s/n^y, r/n^> 

x, and (1/n) Y^l=i ~* Jo ^V- Explicitly evaluating g rs in ( |36*l l 
and taking the n — > oo limit, we obtain the integral eigenvalue 
equation 



dyK(x,y)^ {a \y) = A (q) ^ {a) (x) 



(43) 



with kernel 

K(x,y) = S(x-y) 



1 



7r(a; - y) 



sin [ttA(x - y)] . (44) 



Since the determinant of a matrix is the product of its eigen- 
values, we may write 



E(A) = f[ A (q) = det K. 



(45) 



This is the representation of E as a Fredholm determinant, a 
virtue of which is that n — > oo has been explicitly taken. 

Evaluation of the eigenvalues A^"' is facilitated by a re- 
markable connection between the integral kernel K and 
the prolate spheroidal functions of classical mathematical 
physics. 15 It is found that the prolate spheroidal differential 
operator L and the kernel K commute and therefore have the 
same eigenfunctions. In practice it is easier to determine the 
eigenfunctions of L and then to compute A^ Q ^ by applying K 
to these eigenfunctions. 

Finally, we note that a third expression for E(A) was de- 
rived by Jimbo et al., expressing it as the solution to a second- 
order non-linear differential equation. 30 

Each of these representations has proved useful in the past. 
The asymptotic behavior of E(A) for large A [and hence of 
Pi (A) via Eq. d32b l was derived by des Cloizeaux and Mehta 
by asymptotic analysis of the differential operator L and inde- 
pendently by Dyson by means of an ingenious application of 
inverse scattering theory to the analysis of the kernel K J£ It 
can also be obtained from the Toeplitz representation using an 
asymptotic formula due to Widom.i 5 ^ 

The first numerical computation of E(A) was based on 
evaluation of the Fredholm determinant, taking advantage of 
the connection to prolate spheroidal functionsji 5 - Numerical 
solution of the non-linear differential equation was subse- 
quently found to be more efficient^ 

Here we find that direct evaluation of the Toeplitz determi- 
nant Eq. d40l > for large (but finite) n also provides an accurate 
way to calculate E(A). Fig. [7] shows a plot of E(A) calcu- 
lated in this way with n — 200; the inset shows the conver- 
gence of E(A) as a function of n for A = 1. The results for 
n = 100 and n — 400 differ only in the sixth significant fig- 
ure. The convergence is slower for larger A, but it is clear that 
for all values of A for which the distribution has significant 
weight, and to the accuracy needed for applications to quan- 
tum dots, Toeplitz determinants with n equal to a few hundred 
should suffice. Fig.|7]is the main new result of this subsection. 



B. Consecutive Level Spacing 

In this sub-section we express the consecutive level spac- 
ing (sometimes simply called the level spacing) as a Toeplitz 
determinant. We then show that this expression can be used 
to compute Pi (t) with good precision with modest numerical 
effort. 

For simplicity let us suppose that there are n+2 angles on 
the circle. Pi (t) is the probability that one of these equals t/2, 
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FIG. 7: The spacing determinant E(A) calculated using the Toeplitz 
representation Eq. d40t with n = 200. The inset shows the depen- 
dence of E(A) on n for A = 1. Note the good precision obtained 
with modest numerical effort. 
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FIG. 8: (Color online) The level spacing Pi (A) calculated using 
the Toeplitz representation Eq. J50l > with n = 200 (black points). 
Also plotted is the Wigner surmise (blue curve), known to give an 
excellent approximation to the true distribution. 



another equals —t/2 and the others all lie outside the range 
— 1/2<9 <t/2. Thus we must consider 



P(6i, . . . , 9 n , n+ i — > — , 6 n+2 



(46) 



1 1 



(n + 2)! (2tt) 2 (2tt)' 



it/2 _ —it/2\1 



n 



e it/2 |2| e ie r e -it/2|2 



n 

s—r+l 



Comparing the two equivalent formulations of the circular 
unitary ensemble, Eqs. ( f33l > and Q51 >. we may write 



P(6i, ...,0 n , 9 n +i — * - , #n+2 — * — ^) 



(47) 



1 + 

n 

n 4 [i-c 



9 j 2' 



i - cos {e 3 + -) 



To compute Pi (t) we must now integrate over 9\ , . . . , 9 n out- 
side the range —t/2 < < t/2. Comparing Eq. (07]) to d35l ). 
it is clear that if we define 



2n~t/2 



(16 



6/2 



1 - cos (0 - 1) 



l-cos(0 + -)Jp P (0K(0) (48) 



then 



Pi(t)=/C 1 -^(l-cost)det 5 W. (49) 

The constant /Ciisl/i?i = (27r/n) [see the discussion follow- 
ing Eq. (l32lil. Further, let us work with A, the spacing in units 



of the mean level spacing, rather than t; then 

Pi (A) = lim -4 fl - cos — ") det . (50) 

n^oo V 71/ 

Note that the distribution is normalized and the energy scaled 
so that the mean level-spacing is unity. Eq. ( f50b is the main 
result of this subsection. It expresses the consecutive level 
spacing as a Toeplitz determinant. 

Fig. [8] shows a plot of Pi (A) computed by evaluation of 
the Toeplitz determinant for n = 200 (black points). For com- 
parison we have also plotted the Wigner surmise (blue curve) 
which is known to give an excellent approximation to the true 
spacing distribution. Fig. [8] shows that the n = 200 approxi- 
mation is adequate for any application to quantum dots. 



C. Nearest Neighbor Spacing 

A given level has two neighbors: one above, the other be- 
low. Recently, Forrester and Odlyzkc 25 analyzed the distribu- 
tion of distance to the nearer of the two neighbors, pointing 
out that this distribution is distinct from the consecutive level 
spacing distribution, the traditional object of study in random 
matrix theory. Forrester and Odlyzko adapted the differential 
equation method 30 in their study; here, we show that the dis- 
tribution may also be expressed as a Toeplitz determinant. 

As a preliminary, we introduce a more primitive quantity, 

E nn = n dx 2 ... dx n P(x 1 — > 0,x 2 , ■ ■ ■ ,x n ), (51) 

J out J out 



the probability that one level is at zero energy while the others 
lie outside a band of width 2t. Here 



dx 



dx 



dx 



(52) 



Evidently, i? nn (oo) = and P nn (0) = Pi(0), where P^i is the 
one point correlation function defined in Eq. OTb . 
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By straightforward differentiation we find 

dE nD 



dt 



n(n — 1) / dx$ . 



dx n 



x |p(0, t,x 3 ,..., x n ) + P(0, -t, x 3 ,..., x n )\ 

= i?l(0)P nn (i) (53) 

where the normalization constant ensures that 



poo 

/ dtp nn (t) = l. 

JO 



(54) 



As an application of these identities, consider the Poisson 
model: n independent levels, each distributed uniformly over 
the interval [— n/2, n/2]. In this case it is trivial to calculate 
E nn (t) and, upon differentiation, to find 

lim P nn (£) poisson = 2exp(-2t) . (55) 

n — >oo 

By contrast, the consecutive level spacing is 33 

-PiWpoisson = exp(-i) . (56) 

Thus the Poisson model, which is believed to describe the 
spectral statistics of integrable models, provides an explicit il- 
lustration of the distinction between consecutive level spacing 
and nearest-neighbor level spacing distributions. 

Now let us derive an expression for E nn (t), Eq. ( f5Tb , within 
random matrix theory. Once again, it is convenient to work 
with the circular ensemble, and to suppose that of (n + 1) 
angles, one lies at zero, while the others lie outside the arc 
— t < 9 < t. To this end, we must consider 



(n + i)p(9i,...,e n ,e n+1 

= li ?_fT| e *r 

n\ 2tt (2ir) n 11 1 



0) 



ir n ^ 

l<r<s<n 



(57) 

J9, 1 2 



Comparing the two equivalent forms of the circular unitary 
ensemble, Eqs. d33l and d35l l. we may write 



(n + l)P(0i,. 



,,Q) 



1 1 



J2(-1) P (-1) Q (58) 



jj2^ W (fl i )^ (j) (e i )(l-cose i 



Now, following the reasoning of the two preceding subsec- 
tions leads to the result 



£nn(A) = lim -!-detp e 

n— >oo Z7T 



(59) 



where <? enn is an n x n matrix of the standard Toeplitz form, 
Eq. d4Tb . with the generating function f(8) given by 



fenn(0) 



2(1 - cos 6), for 2^ < 9 < 2tt - ^ 
0, otherwise 



(60) 



The analysis of P nn is entirely similar. The end result is 



P nn (A) = lim — 

n — >oo TL 



1 — COS ■ 



2vrA 



det 



5pnn 



(61) 



where <7 pnn is an n x n Toeplitz matrix generated by 



/cnn(#) 



4(1 



0. 



COS! 



:i-cos[0-^]), 



for 



2ttA 



< e < 2tt 



2ttA 



otherwise 



(62) 

Eqs. (|59l-d62l express E nn (A) and P nn (A) in terms of 
Toeplitz determinants and constitute the main result of this 
subsection. Fig.|5]shows a plot of P nn (A) calculated by eval- 
uation of Toeplitz determinants of size n = 200 (black points). 
The result is in agreement with the earlier calculation by For- 
rester and Odlyzkc 25 based on the method of Jimbo et al. 30 



D. A Pair of Consecutive Spacings 

The joint probability density for two consecutive spacings 
is defined as 

P2(t,£) = Kin{n — l)(n — 2) / dx±... I dx n 

J out J out 

X P{X\ — > 0, X 2 — > T, X 3 — > — I4,...,l n ) 

(63) 

with the excluded region given by 

-e 



dx = / dx + I dx 

out •/— oo 



The normalization constant /Co is needed to ensure that 



/•oo />oo 

/ dr / ^P 2 (r,e) = l. 
Jo Jo 



(64) 



(65) 



To determine the normalization it is convenient to introduce 

E2(j,£,)=nl dx 2 .-- dx n P{xi — > 0, x 2 , ■ ■ ■ , x n ). 

J out J out 

(66) 

Evidently, P 2 (oo, f) = P 2 (r, oo) = and P 2 (0, 0) = Pi(0), 
the one point correlation. It is also easy to see that 



(67) 



and hence /C 2 = l/Pi(0). 

Again, in practice it is more convenient to work with the 
circular unitary ensemble and to compute the probability that 
of (n + 1) angles, one lies at 9 = 0, another at r, a third at — £ 
and none of the others lie on the arc — £ < 9 < r. It is also 
convenient to define r = 27rA 2 /n and £ = 27rA 1 /r7,; Ai and 
A 2 are the level spacings in units where the mean spacing is 
one. 
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FIG. 9: (Color online) The joint spacing distribution P2(Ai, A2) in 
the GUE, calculated using the Toeplitz representation Eq. d68t with 
n = 200, as a function of Ai while A2 is held fixed at 0.5 (red 
points), 1.0 (green points), and 2.0 (blue points). The curves corre- 
spond to the generalized Wigner surmise approximation, Eq. (16). 



The calculations closely parallel those in the preceding sub- 
sections. The end result is that 



P 2 (A 1; A 2 ) = lim — 

n—»oc fi° 

'2ttA 2 



1 — cos 



(68) 



1 _ cos (^)] [l - cos { — (A x + A 2 )}] det 3 ( 2 > 



Here g( 2 ' is an n x n Toeplitz matrix generated by 



cos I 



/cnn(#) 



8[l-cos6>] [l 

x[l-cos(0-f A 2 )], 
for ^A 2 < 9 < 2tt 
0, otherwise 



: A X 



(69) 

Fig.|9]is aplot of P 2 (Ai, A 2 ) calculated using Eq. ( 168) with 
n = 200. Fig.|9]and the Toeplitz representation of P 2 (Ai , A 2 ) 
in (l68l l and ( |69l are the main results of this subsection. Mehta 
gives an expression for P 2 in terms of eigenvalues and eigen- 
functions of the kernel, Eq. (|44) . and tabulates P 2 for a small 
number of (Ai, A 2 ) values^! Our numerical scheme allows 
us to replicate and extend those results with only modest com- 
putational effort. 



E. Asymptotics 

In this sub-section we employ the Toeplitz representation 
to obtain large A asymptotic expressions for P nn (A) and 
P nn (A), analogous to the known results 15 for E(A) and the 
consecutive level spacing Pi (A). From a practical point 
of view these expressions are not very useful because the 
level spacing distributions have negligible weight in the tails. 
Nonetheless they are of some theoretical interest as they rep- 
resent one limit in which it is possible to obtain exact analytic 
results for the level spacing problem. Moreover, it is precisely 



when A ^> 1 that numerical computation of the distribution 
is most difficult. 

The asymptotic results are obtained by straightforward ap- 
plication of a formula due to Widom. 32 Consider an n x n 
Toeplitz matrix T with generating function f(9) non-zero and 
positive on its support a < 9 < 2tt — a, and zero outside this 
arc. According to Widom, for large n 

detT w [nsin-j [cos -J (70) 
x exp(26) exp(2P ) exp f i kF k F^ k ^] . 

Here 



k=l 



b = — ln2+ -C'(-l) ~ -0.219250583, (71) 
24 2 



and the function 

Fie) = 



(0) = f ( 



2 cos 



a 

cos — cos ( 
2 



(72) 



is defined over the the range < 9 < 2tt. Fq and F k are Fourier 
coefficients in the expansion 



hxF(6) = F k exp(ik9). 



(73) 



k— — c 



For our application f(8) = 2(1 — cost?) and a = 27rA/n 
with 1< A < n [Eqs. (g9) and d60"Vl, Thus 



(74) 



F(9) =4(1- cos^ - cos^ 9 



The Fourier coefficients of In F may be determined by con- 
tour integration around the unit circle. We find 



Fn = 4 In ( cos - J - 2 In \1 - sin - 



1 - sin(a/2) 
cos(a/2) 



|fc| 



(75) 



The expression for F k applies provided k is non-zero and 
even. For odd k, F k = 0. 

Substituting Eq. (fTBT l in (l70b . performing the fc-sum, noting 
that a = 2ir A/n, and taking the large n limit, we finally obtain 



- det g enn 
n 



(76) 



exp 



7T 2 5 

A 2 + 2vrA - - InOA) + 2b - In 4 

2 4 



E nn follows directly from Eq. d59l ), and Eq. ((53) then implies 
that in the large n limit 



Pnn(A) 

from which we obtain 



dA 



P nn « I ex P (26)^ T7I 



exp 



• det g c 



ir 2 A 2 



2vrA 



(77) 



(78) 



13 



Eqs. d76l l and (1781 are the asymptotic formulae that we 
sought. We see that P lm (A) has a much more rapid decay for 
large A than the consecutive level spacing; the latter is known 
to have the leading behavior Pi (A) oc exp(— 7r 2 A 2 /8). This 
more rapid decay is already evident at quite modest values of 
A upon inspection of Figs. |5]and[8] 

V. SUMMARY AND CONCLUSION 



suits was to test the simple approximate formulae developed 
in Section [TTTJ but the exact results may be of intrinsic inter- 
est from the viewpoint of fundamental random matrix theory 
also. 

The derivation of the universal Hamiltonian model as- 
sumes that the wavefunctions are distributed according to a 
pure ensemble of random matrix theory. The problem of 
electron-electron interactions in the crossover between sym- 
metry classe s 34 ' 35 ' 36 ' 37 ' 38 ' 39,40 ' 41 needs careful reconsideration. 



In this paper we have developed the relation between ran- 
dom matrix theory and electron-electron interaction in two di- 
mensional electron gas quantum dots within the framework 
of the universal Hamiltonian model. We find the probability 
distribution of the spacing between conductance peaks in the 
Coulomb blockade regime for the standard orthogonal or uni- 
tary cases (the symplectic case is trivial, being essentially the 
same as the constant interaction model). For these cases there 
are departures from the constant interaction model even in the 
paramagnetic limit where J<C (A). In this limit, the Coulomb 
blockade peak spacing distribution can be expressed in terms 
of the joint probability distribution of a few consecutive level 
spacings within RMT Making use of approximate formulae 
for the random matrix joint probability distributions, we ob- 
tain analytic formulae for the Coulomb blockade peak spacing 
distribution. These expressions reproduce earlier numerics 20 
and should facilitate future comparisons to experiments done 
at low temperature (k B T < 0.1(A)) 

The simple approximate formulae for the joint probability 
distributions P^Ai, A2) and p3(Ai, A2, A3) are obtained in 
Section [HI] by analogy to Wigner's famous surmise for the 
consecutive spacing distribution Pi (A). In Section HVl we 
have also derived a number of new exact results in random ma- 
trix theory. These include a representation of the consecutive 
level spacing distribution Pi (A), the nearest neighbor spac- 
ing P nn (A), and the joint spacing distribution p2(Ai, A2) 
as Toeplitz determinants. The Toeplitz representation enables 
efficient numerical computation of these distributions and al- 
lows us to infer the previously unknown large A asymptotics 
of P nn (A). Our primary motivation in deriving these exact re- 
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APPENDIX A: PEAK SPACING DISTRIBUTION 
INCLUDING S= § 

The contribution to the probability distribution of the 
Coulomb blockade peak spacing of each spin transition is 
obtained by using the following standard procedure: if s — 
f(x, y, z) is a function of random variables x, y, and z, then 
its probability distribution is given by 



P(s)- 



P(x,y,z)5(s-f(x,y,z))dxdydz (Al) 



where P(x,y,x) is the joint probability distribution of the 
three random variables. Now, since in our case the function 
f(x, y, z) has different forms for different values of the vari- 
ables (see Table we analyze each case separately and add 
the appropriate step functions (9). The total spacing distri- 
bution is then the sum of each contribution. In what follows, 
both the even and the odd distribution are normalized to one. 



1. Even distribution 



\x M) {s) = @(s-\j) J J P 3 (x,a+p ) z)e(8+x~j)Q(a+ 

p ( i,i,i ) ( s )=e( s --j)e(-j- s ) J J p 3 (x,-j- s ,z)dxdz 



Z— — J) dxdz 



(A2) 
(A3) 



1 ,00 ,min(2J,§J-s) 1 \ / \ 

P (i,i,|)(s) = 0(s+2 J J i / o dx J ft(x,y,8+-J)Q[x+y-3J)dy 



(A4) 



j.3J pmin(2J,3J-z) g 1 \ / 5 \ 

P ( | )liS) (s) = y o dzj P 3 {-J+s-y-z,y,z)e{z-s+-J)Q(s+-J-y-z)dy (A5) 
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2. Odd distribution 



/ 3 \ f°° f°° 

P(o^,o)(s) = S[s--J) / / P 2 (x,y)dxdy (A6) 
P(i,i,o)( s ) = @ ( s +l J ) @ (l J - s ) I P*(s+\j,y)e(s+y-h) dy (A7) 

«* 2 J 

P{i,\, i)(*)=e(* ~5 J )/ ft(:M+^-z)e(s-»-~j)e(|j+a-z)da: (as) 

^(i,f,i)( s ) = e ( s -5 J ) / ft(aJ,^-»-a:)e^J- a -x)e(a;+8-|j)dx (A9) 

3. Explicit expressions for the GOE case 
In the GOE case (see Fig. [0, the integrals above can be done explicitly. For the even distribution we get 

P (f°o,i)( s ) = f Q { s ~\ J ) j©(s-^)4xe- 3a;2 (Xx-e* 2 (2x 2 ~l) Erfc(x)) +e(~ J-a^e" 3 ^ 2 ^-^ 

^(2y + x)(l + 2y{y+x)) + e {2y+x)2 (l + 2y 2 {l+y 2 ) + 2yx (l + 2y 2 )-2x 2 (l + y 2 ) -Ayx 3 ) Erfc(2y+.x) \ J 

(A10) 

wherea;=(s+|j)//V3and2/=(|j-s)/7v / 3, and 

p ( f, e ii)( s ) = /e( s -ij)e(^j- s )4^- 3 ^ 2+2 ^ 2 )^(2y+.x)(i+2j/(y+x)) 

+e (2y +a; ) 2 (i + 2y 2 (l+t/ 2 ) + 2j/a;(l + 2j/ 2 )-2a; 2 (l + y 2 ) -4ya; 3 ) Erfc(2y+x)^j (All) 

with x = (| J — s)//\/3 and y = (i J+s)f The corresponding expressions for Eqs. (lA4b and (IA5l l are similar in form but 
too cumbersome to be presented here. 

The partial contributions to the odd distribution are given by 

P^ ifi )( s )= S { S -l J y j2 (^" 3j2 "(2j 2 -l)Erfc(3^)^ (A12) 

with j = J\f2h, and 

P (T1 o)( s ) = VW (j+x)e ~ 4( ^ )2 {' ' (nnj~x)-e^^ 2 7ri (2( J+ x) 2 -l)Erfc(ll^x)) e(^J-s) 



x6(ij+s) +e~ 16 ^ 5 ^ (2(9j+x)-e^ j ^ 2 ^ (2(j + x) 2 -l) Erfc(9j+z)) &{s~ j)e(|j-s) 



(A13) 



P ( ^ ii)(s ) = ^(5^x)e- 4 ( 21 ^"^ 2 ) (2(3 J - 2 ;) + e ^ 2 7r^ (2(5j + x) 2 -l) Erf(3j-x)) ( S -1 j)e(|j- S ) 



(A14) 
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P ( s° 3 1} ( S ) = ^l^(7j-x)e- 4 ^ 2 ~ w ^^ je- 48 -'' 2 (2(j+x)+e^\i (2(x-7j) 2 -l) Erf(j + a;)) © 

x6(Jj- S ) +e 16 ^ (2(7j-x) + e ^) 2 7r5 (2(:r-7j) 2 -l) Erf(7j-a;)) e( s -|j)e(^ J- S )| 



(A15) 



with j = J y / b/2A and x = s 
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